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Abstract. We consider a class of asset pricing models, where the risk-neutral 
joint pr ocess of log-price and its stochastic variance is an affine process in the 
sense of iDuffie. Filipovic. and Schachermaverl [20031 ]. First we obtain condi- 
tions for the price process to be conservative and a martingale. Then we 
present some results on the long-term behavior of the model, including an ex- 
pression for the invariant distribution of the stochastic variance process. We 
study moment explosions of the price process, and provide explicit expressions 
for the time at which a moment of given order becomes infinite. We discuss 
applications of these results, in particular to the asymptotics of the implied 
volatility smile, and conclude with some calculations for the Heston model, a 
model of Bates and the Barndorff-Nielsen-Shephard model. 



1. Introduction 



Duffie. Pan, and Singletonl [2000( | introduced the notion of an affine jump-diffusion, 



which is a jump-diffusion process, whose drift vector, instantaneous covariance ma- 
trix and arrival rate of ju mps all depend in an affine way on the state vector. 
Duffie. Pan, and Singleton! remark that models built on affine processes provide a 
balanced tradeoff between analytical tractability and complexity, making them an 
attractive choice for applications in mathematical finance. In particular they men- 
tion applications t o the pr i cing of options in stochas tic vo latility models and note 
that the models of lHestonl [l993| . iBatesl [l996l [2000j ] , and lBakshi. Cao. and Chen 
19971 fall into the affine class. To these, we could add the more recent models of 



Barndorff-Nielsen and Shephardl [2001 1 and Carr and W u 2004] , which are also of 
affine type. 



Duffie. Filipovic. and Schachermaverl [20031 ] subsequently extended the class of affine 
jump-diffusions, defining an affine process as a time-homogenous Markov process, 
whose characteristic function is the exponential of an affine function of the state 
vector. It turns out that this class coincides for a large part with the class of affine 
jump-diffusions, but al so allows for infinite activity of jumps and for killing or ex- 
plosions of the process. IDuffie. Filipovic. and Schachermaverl aim to give a rigorous 
mathematical foundation to the theory of affine processes, covering many aspects, 
such as the characterization of an affine process in terms of the 'admissible param- 
eters' (comparable to the characteristic triplet of a Levy process) and properties of 
the ordinary differential equations ('generalized Riccati equations') that are implied 
by the process. 



Key words and phrases, affine process, stochastic volatility, moment explosions, implied volatility 
smile. 

Supported by the Austrian Science Fund (FWF) through the START programm Y328. 



2 



MARTIN KELLER-RESSEL 



In this article we study stochastic volatility models, comprised of a log-price 
process {Xt)t>o and a stochastic variance process (Vt)t>o, such that the joint pro- 
cess (X t , Vt)t>o is an affine process. We will show that many properties of such a 
model, including its long-term behavior and moment explosions, can be analyzed 
by studying differential equations of the generalized Riccati type. Our results on 
the long-term behavior are formulated as asymptotic results for the cumulant gen- 
erating function of the stock price, a s time goes to infinity. Asymptotics of this 
type have been used by Lewis! 2000| | to obtain large-time-to-maturity results for 



the implied volatility smile of stochastic volatility models via a saddlepoint expan- 
sion. The issue of moment explosions in stoc hastic volatility models has rec ently 
received much attention, due to the articles of Andersen and Piterbar and 



Lions and Musiela |2007l |. Moment explosions are intimatel y connecte d to large- 



strike asymptotics of th e implied volatility smile via results of [Leg [2004j , that have 
later been expanded by Benaim and Frizf 2006j . 



In the first part of the paper we introduce our main assumpt ion, that the .join t 



process (Xt, Vt)t>o is affine, and recapitulate the main results of lDuffie et al.l [20031 ] . 
We derive necessary and sufficient conditions for conservativeness of the process 
and for the martingale property of the discounted price process St = exp(A t ). At 
the end of Section [2] we add two assumptions, and give a precise definition of the 
class of affine stochastic volatility models, which constitutes the main subject of this 
article. In Section[3]we derive our central results on long-term properties of an affine 
stochastic volatility model, providing conditions for the existence of an invariant 
distribution of the stochastic variance process, and characterizing this distribution 
in terms of its cumulant generating function. We also give results on the long-term 
properties of the price process, showing that as time tends to infinity, the marginal 
distributions of the price process approach those of an exponential-Levy process. 
The characteristic exponent of this Levy process can be derived directly from the 
specification of the affine stochastic volatility model. Both results are obtained by 
applying qualitative ODE theory to the generalized Riccati equations introduced 
in the first part. 

In Section [4] we study moment explosions of the price process, and show that an 
explicit representation for the time of moment explosion can be given - not only 
for the primary model, but also for the model in the stationary variance regime. 
In Section [5] we outline applications of our results to the asymptotics of im plied 
volati lities and of implied forward volatilities. We briefly discuss the results of iLeel 
[2004j and point out the connection between the stationary variance regime and 
the pricing of forward-start options, when the time until the start of the contract 
is large. We conclude in Section [6] with explicit calculations for several models to 
which our results apply, such as the Heston model, a Heston model with added 
jumps, a model of iBatesl . and the Barndorff-Nielsen-Shephard model. 



2. Affine Stochastic Volatility Models 

2.1. Definition and the generalized Riccati equations. We consider an asset- 
pricing model of the following kind: The interest rate r is non-negative and constant, 
and the asset price (St)t>o is given by 



S t = exp(ri + A t ) t>0, 
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such that (Xt)t>o is the discounted log-price process starting at Xq 6 M a.s. The 
discounted price process is simply exp(JTt), such that we will assume in the remain- 
der that r = 0, and that (St)t>o is already discounted. Denote by (Vt)t>o another 
process, starting at Vq > a.s., which can be interpreted as stochastic variance 
process of (X t ) t >o, but may also control the arrival rate of jumps. The following 
assumptions are made on the joint process (X t , Vt)t>o- 

Al: (Xt, Vt)t>o is a stochastically continuous, time-homogeneous Markov pro- 
cess. 

A2: The cumulant generating function $t(u, w) of (X t , Vt) is of a particular 
affine form: We assume that there exist functions (f)(t,u,w) and ip(t,u,w) 
such that 

$t(u,w) := logE[exp(uX t + wV t )\X ,V ] = cj)(t,u,w) + V ip(t,u,w) + X u 

for all (t,u,w) € M>o x C 2 , where the expectation exists. 
By convention, the logarithm above denotes the principal branch of the complex 
lo garithm. Assumpti ons Al and A2 make (X t , Vt)t>o an affine process in the sense 
of lDuffie et all [2003j . The term Xqu m the cumulant generating function <fr t (u,w) 



corresponds to a reasonable homogeneity assumption on the model: If the starting 
value Xq of the price process is shifted by x, also X t is simply shifted by x for any 
t > 0. Note that Assumption A2 also implies that the variance process (Vt)t>o is a 
Markov process in its own right. We do not yet make the assumption that (St)t>o 
is conservative (i.e. without explosions or killing) or even a martingale. Instead it 
will be our first goal in Section 12.21 to obtain necessary and sufficient conditions for 
these properties. 

Applying the law of iterated ex pectations to $>t(u ,w) yields the following 'flow- 
equations' for <j> and ip: (see also lDuffie et all [20031 Eq. (3.8)-(3.9)]) 

<j>(t + s, u, w) = 4>{t, u, w) + <j)(s, u, ip(t, u, w)), 

^ ^ ip(t + s,u,<w) = ip(s,u,ip(t,u,w)), 

for all t,s > 0. The following result will be crucial: 

Theorem 2.1. Suppose that \<j)(T,u,j])\ < oo and \ip[r,u,rj)\ < oo for some 
(t, u, rj) € M^o x C 2 . Then, for all t £ [0, r] and w £ C with Re w < Re rj 

\<f)(t, u, w)\ < oo, \<f)(t, u, w)\ < oo, 

and the derivatives 
(2.2) F(u,w) := — <j)(t,u,w) 



t=o+ 



d 

R(u,w) := —ip(t,u,w) 



t=o+ 



exist. Moreover, for t £ [0, t), <f> and ip satisfy the generalized Riccati equations 
(2.3a) d t (f)(t,u,w) = F(u,ip(t,u,w)), (f>(0, u, w) = 

(2.3b) dtip(t,u,w) = R(u,ip(t,u,w)), ip(0,u,w)=w. 

The above theorem is 'essentially' proven in lDuffie et alj [20031 ] . but under slightly 
different condition^. Note that the differential equations (|2.3[) follow immediately 
from the flow equations (12. ip by taking the derivative with respect to s, and evalu- 
ating at s — 0. They are called generalized Riccati Equations since they degenerate 



] Duffic et alj assume differentiability of cj> and tp with respect to t ('regularity') a priori, while in 
our case we can deduce it directly from Assumption A2. A proof is given in the appendix. 
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into (classical) Riccati equations with quadratic functions F and R, if (Xt,Vt)t>o 
is a pure diffusion process. 

Note that the first Riccati equation is just an integral in disguise, and 4> may be 
written explicitly as 

(2.4) (j)(t,u,w)= / F(u,tp(s,u,w)) ds . 



Also the solution ip of the second Riccati equation can be represented at least 
implicitly in the following way: Suppose that tp(t,u,w) is a non-stationary local 
solution on [0, S) of (|2.3bj) . Then R(u, ip(t, u, w)) ^ for all t e [0, 5), and ip(t, u, w) 
is a strictly monotone function of t; dividing both sides of (|2.3b[) by R(u, tp(t, u, u>)), 
integrating from to t < 5, and substituting r/ = ip(s, u, w) yields 

ri>(t,u,w) ^ 



(2.5) / 



-ds = t 



R(U,T)) 

Another important result that can be found in lDuffie et al. .2003] states that F 
and R must be of Levy-Khintchine form, i.e. 

(2.6a) F(u,w) = (u,w)-^-(l\+b-(l 



D\{0} 



e^-l-^f^J-Ml m(dx,dy) 



(2.6b) R(u,w) = (u,w)-^-r\+P-r) , 



D\{0} 



1 - uj R (x,y) ■ ( " I I ii(dx,dy) 



where D = lx R^Oi an( i ^f, ^>r are suitable truncation functions, which we fix 
by defining 







uJ F (x,y) = Hf and u R (x,y) = 1+ y x 



1+3/2 



Moreover the parameters (a, a, b, (3, c, 7, m, [i) satisfy the following admissibility 
conditions: 

• a, a are positive semi-definite 2 x 2-matrices, and an = 0,21 = &22 = 0. 

• b € D and (3 e R 2 . 

• 0,7 S R^o 

• m and /i are Levy measures on D, and /d\{o} (( a;2 + f ) ^ m {dx, dy) < 
00. 

The affine form of the cumulant generating function, the generalized Riccati 
equations and finally the Levy-Khintchinc decomposition (|2.6[) lead to the follow- 
ing interpretation of F and R: F characterizes the state-independent dynamic of 
the process (X t ,V t ) while R characterizes its state-dependent dynamic. Both F 
and R decompose into a diffusion part, a drift part, a jump part and an instan- 
taneous killing rate. Hence a + aV t can be regarded as instantaneous covariance 
matrix of (X t ,V t )t>o, b+ V t (3 as the instantaneous drift, m(dx,dy) + V t fi(dx,dy) 
as instantaneous arrival rate of jumps with jump heights in (dx x dy), and finally 
c + "fVt as the instantaneous killing rate. 
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The following Lemma establishes some important properties of F and R as func- 
tions of real-valued arguments. A proof is given in the appendix. 

Lemma 2.2. (a) F and R are proper closed convex functions on M 2 . 

(b) F and R are analytic in the interior of their effective domain. 

(c) Let U be a one- dimensional afftne subspace ofM 2 . Then F\jj is either a strictly 
convex or an affine function. The same holds for R\u ■ 

(d) If (u, w) G domF, then also (u, 77) € domF for all n < w. The same holds for 
R. 

Remark 2.3. As usual in convex analysis, we regard F and R as functions defined 
on all of R 2 , that may attain values in lU {+00}. The set {(u, w) : F(u, w) < 00} 
is called effective domain of F, and denoted by domF. 



We define a function x(it), that will appear in several conditions throughout this 
article. Corollary 13.51 gives an interpretation of x as a rate of convergence for the 
asymptotic behavior of the cumulant generating function of (X t )t>o- 

Definition 2.4. For each u £ R where R(u, 0) < 00, define x{ u ) as 

X(u) := ^(u,w) 

x(u) is well-defined at least as a limit as w | 0, possibly taking the value +00; 
it can be written explicitly as 

f 

/d\{o> " V 1 + y' 2 

Note that also x( u ) 1S a convex function. 



X(u) = a 12 u + Pi + I y[e xu -— — -)n{dx,dy) 



2.2. Explosions and the martingale property. We are interested in conditions 
under which St = exp(X t ) is conservative and a martingale. If such conditions are 
satisfied, (S t )t>o may serve as the price process under the risk-neutral measure in 
an arbitrage-free asset pricing model. The following theorem gives sufficient and 
necessary conditions: 

Theorem 2.5. Suppose (X t ,Vt) satisfies Assumptions Al and A2. Then the fol- 
lowing holds: 

(a) (S t )t>o is conservative if and only if F(0,0) = i?(0,0) = and 

dn 



(27) Jo-mv) 

(b) (S t )t>o * s a martingale if and only if it is conservative, F(1,0) = i?(l,0) = 
and 

Remark 2.6. The notation J Q _ denotes an integral over an arbitrarily small left 
neighborhood of 0. 

By flUt the condition F(0, 0) = #(0,0) = is equivalent to c = 7 = 0, i.e. 
obviously the killing rate has to be zero for the process to be conservative. As 
will be seen in the proof, the integral conditions (|2.7p and l[2.8[l are related to a 
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uniqu eness condition for non-Lipschitz ODEs, which has been discovered by lOsgoodI 
[l898j . 

The following Corollary gives easy-to-check sufficient conditions: 
Corollary 2.7. Suppose (Xt,Vt) satisfies Assumptions Al and A2. 

(a) If F(0, 0) = R(0, 0) = and x(0) < oo then (St)t>o * s conservative. 

(b) If (St)t>o is conservative, F(1,0) = i?(l,0) = and %(1) < oo, then (St)t>o is 
a martingale. 

Proof. For a proof of [2~3a1 we refer to [Filipovid . [200l Th. 4.11]. Statement l2ilEl 



can be shown in a similar way: 

Since (X t , Vt) is Markovian, we have for all < s < t, that 

E[S t \F a ] - S,exp(0(t-a > l,O) + V , .V(t-*,l,O)) . 
We have assumed that Vb > a.s., such that (St)t>o is a martingale if and only if 
(Xt)t>o is conservative and ip(t, 1, 0) = 0(t, 1, 0) = for all t £ K^o- 
We show Corollary 12.71 and the first implication of 12. jbl Suppose that (St)t>o is 
conservative and that F(1,0) = R(1,Q) = 0. By Theorem 12. II ibit. l,w) solves the 
differential equation 

(2.9) ^(t,l,w)=R(l,^(t,l,w)), ^(0,l,w)=w 

for all w < 0. Since R(l, 0) = it is clear that ^(i, 1,0) = satisfies this ODE for 
the initial value w = 0. To deduce that ^(t, 1,0) = ip(t, 1, 0) however, we need to 
know whether the solution is unique. Since R(l,w) is continuously differentiable 
for w < 0, it satisfies a Lipschitz condition on (— oo, 0). If %(1) < oo, the Lipschitz 
condition can be extended to (— oo,0], and ip(t, 1,0) = is the unique solution. 
Without the assumption that x(l) < oo, we substitute Lipschitz' condition by 
Osgood's conditionj (|2.8|) : Suppose that (|2.8|) holds, and there exists a non-zero 
solution ip such that ift(ti, 1, 0) < for some t\ > 0. Then for all t < t\ such that 
ip remains non-zero on [t,ti] we have (similarly to (|2.5p ) that 

(2.10) / -P— =t-t 1 . 



V>(ti,i,o) 



R(i,n) 



Assume that to > is the first point left of t\ such that t/)(to,l,w) = 0. Letting 
t i to, the left side of (|2.10p tends to — oo, whereas the right side remains bounded, 
leading to a contradiction. We conclude that ip(t, 1, 0) = is the unique solution of 
(123)) . Finally equation flggp together with F(l, 0) = yields that also <j)(t, 1, 0) = 
for all t £ M^o and we have shown that (St)t>o is a martingale. 

For the other direction of I2.db1 note that (St)t>o being a martingale implies 
that (f> = i/j = solve the generalized Riccati equations and thus that F(1,0) = 
R(1,0) = 0. It remains to show ([231) . Assume that (|2~8]) does not hold. Then, 
for each t > 0, (|2.10p with t\ = implicitly defines a solution ijj(t, 1,0) of the 
generalized Riccati equation (|2.9p . satisfying ^>(£, 1, 0) < for all t > 0. By unique- 
ness of the solution ip(t,l,w) for w < and the flow property (|2.ip . we have 
ijj(t + s, 1,0) = V , (£;>1>V'( S )1>0)) for i, s small enough. Letting s | we obtain 
^(i, 1, 0) = ip(t, 1, 0) < 0, which is a contradiction to ip = 0. □ 



2 See lOseoodl |l898H 
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We add now two assumptions to Al and A2 and complete our definition of an 
affinc stochastic volatility model: 

A3: The discounted price process St — e Xt is a martingale. 
A4: R(u, 0)^0 for some «£i, 

Assumption A4 excludes models where the distribution of (Xi)t>o does not de- 
pend at all on the volatility state Vo- In such a case we can not speak of a true 
stochastic volatility model, and it will be beneficial to avoid these degenerate cases. 
We are now ready to give our definition of an affine stochastic volatility model: 

Definition 2.8. The process (X t ,Vt)t>o is called an affine stochastic volatility 
model, if it satisfies assumptions Al - A4. 

A simple consequence of this definition, that will often be used is the following: 

Lemma 2.9. Let (A t , Vt)t>o be an affine stochastic volatility model. Then R(u,0) 
is a strictly convex function, satisfying R(0, 0) = i?(l,0) = 0. 

Proof. From assumption A3 and Theorem 12.51 it follows that i?(0, 0) = i?(l,0) = 
0. Lemma 12.21 implies that 0) is either strictly convex or an affine function. 
Assume it is affine. Then R(u, 0) = for all «eR. This contradicts A4, such that 
we conclude that R(u, 0) is a strictly convex function. □ 

3. Long-term asymptotics 

In this section we study the behavior of an affine stochastic volatility model as 
t — > oo. We focus first on the stochastic variance process (Vt)t>o- Under mild 
assumptions this process will converge in law to its invariant distribution: 

3.1. Stationarity of the variance process. 

Proposition 3.1. Suppose that Al and A2 hold, that %(0) < and the Levy 
measure m satisfies the logarithmic moment condition 

/ (logy) m(dx,dy) < oo. 

Jy>l 

Then (Vt)t>o converges in law to its unique invariant distribution L, which has the 
cumulant generating function 

iKeller -Ressel and Steinerl |2008| show that under the given conditions the pro- 
cess (Vt)t>o converges in law to a limit distribution L, whose cumulant generating 
function can be represented by ()3. 1|) . A short argument at the end of this para- 
graph shows that the limit distribution is also the unique invariant distribution 
of (Vt)t>o- First we make the following definition: Given some affine stochastic 
volatility model (X t ,Vt)t>o, we introduce the process (X t ,Vt)t>o, defined as the 
Markov process with the same transition probabilities as (Xf,Vt)t>o, but started 
with Xq = and Vq distributed according to L. We will refer to (X t , Vt)t>o as the 
stochastic volatility model (X t , Vf)t>o 'in the stationary variance regime'. We also 
define the associated price process St := exp(rt + X t ). As we discuss in Section [5] 
the process (X t ,Vt)t>o can be related to the pricing of forward-starting options, 
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when the time until the start of the contract is large. 
The cumulant generating function of (X t , Vt) is given by 



(3.2) 
log I 



e 



uXt+wV t ] 



logE 



exp 



w) + Voip(t, u, w) 



4>{t, u, w)+l(ip(t, u, w)) 



We verify now that L is indeed an invariant distribution of (Vt)t>o : 



(3.3) E exp(wVt\ = exp (<f>(t, 0, to) + l(ip(t, 0, w))) 



rt r>< 

exp | / F(0,%p(s,0,w))ds + 



rif>(t,0,w) 



exp 



F(0,y) 



dn + 



o 



V>(t,0,' 



R(Q,rj) 

R(o,v) 



dq 



d-q = exp(Z(to)) 



where we have used that under the conditions of the Propo sition above, t^(t, 0, w) is 
a stric tly monotone function converging to as t — > oo. (cf. iKeller-Ressel and Steiner 
2008j l. To see that L is unique, assume that there exists another invariant dis- 
tribution L', and let (y/) t > be the variance process started with Vq distributed 
acco rding to L' . Again we use that <j>(t, u, w) — * l(w) and ip(t, 0, w) ~+ as t — > oo 



fsee IKeller-Ressel and Steiner! [2008]), and get that 



lim E [exp(wV t ')] = E lim exp (</>(*, 0, w) + V V(t, 0, w)) = E[exp(/(u;))] = e l{w) 
for all w < 0, in contradiction to the invariance of L' . 

3.2. Long-term behavior of the log-price process. We have seen that (Vt)t>o 
converges to a limit distribution, but we do not expect the same for the log-price 
process (JT t ) t > . Nevertheless, it can be shown that the rescaled cumulant gener- 
ating function |logE [e x * n ] converges under suitable conditions to a limit h(u), 
that is again the cumulant generating function of some infinitely divisible random 
variable. This result can be interpreted such, that for large t the marginal distri- 
butions of (Xt)t>o are 'close' to the marginal distributions of a Levy process with 
characteristic exponent h[u). Furthermore, h(u) can be directly obtained from the 
functions F and R, without knowledge of the explicit forms of <j> and ip. We start 
with a preparatory Lemma: 

Lemma 3.2. Let (X t ,Vt) t >o &e a- n °ff^ n ^ stochastic volatility model and suppose 
that x(0) < and x(l) < 0. Then there exist a maximal interval I and a unique 
function w G C(I) fl C 1 (/°), such that 

R(u,w(u)) = for all u G I 

and w(0) = w(l) = 0. 

Moreover it holds that [0,1] C I, w(u) < for all u £ (0,1); w(u) > for all 
u £ I\ [0, 1]; and 

(3.4) 



^(u,w(u)) < 



for all u£ 1° 



We show Lemma I3~2l together with the next result, which makes the connection 
to the qualitative properties of the generalized Riccati equations. 
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Lemma 3.3. (a) For each u G 1° , w(u) is an asymptotically stable equilibrium 
point of the generalized Riccati equation (|2.3bj) . 

(b) For u G 1° , there exists at most one other equilibrium point w{u) ^ w(u), and 
if it exists, it is necessarily unstable and satisfies w(u) > max(0, w(u)). 

(c) For u G M \ /, no equilibrium point exists. 

Proof. Define L — {(u,w) : R(u,w) < 0}. As the level set of the closed con- 
vex function R, it is a closed and convex set. For all u G R, define w(u) — 
inf {w : (u,w) G L}, and / = {u G R : w(u) < oo}. Clearly w(u) is a continuous 
convex function, and / a subinterval of R. We will now show that w(u) and / 
satisfy all properties stated in Lemma 13.21 By assumption A3 and Theorem 12.51 
i?(0, 0) = R(l, 0) = 0; together with Lemmal2T2lit follows that the set [0, 1] x (-oo, 0] 
is contained in domi?. Since R(u,0) is by Lemma 12.91 strictly convex, and also 
x( u ) is convex, we deduce that R(u, 0) < and §^(u 7 0) = xi u ) < f° r au 
u G (0, 1). In addition R(u, w), as a function of w, is either affine or strictly convex, 
such that there exists a unique point w(u), where R(u,w(u)) = 0, and necessarily 
jf^{u,w(u)) < 0. It is clear that for u G (0,1) w(u) coincides with the function 
defined above, and that w(u) < 0. At u = we have that R(0, 0) = and x(0) < 0, 
implying that w(0) = 0. A symmetrical argument at u = 1 shows that w(l) = 0, 
and thus that [0, 1] C I. 

We show next that w(u) G C 1 (/°): Define u + = sup/, and w + = lim„| n+ w(u); 
U—,W— are defined symmetrically at the left boundary of /. Note that u± and w± 
can take infinite values. Define the open set 

K := {(Ait_ + (1 - X)u + , w) : A G (0, 1), w < Aw_ + (1 - X)w+} . 

Lemma 12.21 implies that K is contained in the interior of dom R. On the other 
hand, the graph of w, restricted to 1° . i.e. the set {(u,w(u)) : u G 1°}, is clearly 
contained in K. Since R is by Lemma 12.21 an analytic function in the interior of 
its effective domain, the implicit function theorem implies that w(u) G C 1 (J°). 
In addition it follows that ^(u,w(u)) ^ for all u G 1° , such that the asser- 
tion §^j(w, w(u)) < 0, which we have shown above for u G (0, 1), can be extended 
to all of 1° . The claim that w(u) > for u G I \ [0, 1] can easily be derived from 
the convexity of w(u), and the fact that w(u) < inside (0, 1) and w(0) — w(l) = 0. 

We have now proved most part of Lemma 13.21 (except for the uniqueness) , and 
turn towards Lemma 1331 Since R(u, w(u)) — and §^(u, w(u)) < for all u G 1° , 
w(u) must be an asymptotically stable equilibrium point of the generalized Riccati 
eauation l2.3bl showing |3. 3b . Assume now that for some u G 1° there exists a point 
w(u) ^ w{u) such that R(u,w(u)) = 0. By Lemma [2.21 R(u,w) is, as a function 
of w, either strictly convex or affine. If it is affine, it has a unique root, and w(u) 
cannot exist. If it is strictly convex, there can exist a single point w(u) other than 
w(u), such that R(u,w(u)) — 0. Necessarily w{u) > w(u) and ^(u,w(u)) > 0. 
This shows that w(u) is an unstable equilibrium point of the generalized Riccati 
equation for In addition w(u) > w(u), and in particular the fact that w(0) > 
and w(l) > shows the uniqueness of w(u) in the sense of Lemma 1331 To see that 
w(u) > max(0, w(u)), note that we only have to show that w(u) > 0, whenever 
w(u) < 0. This is the case only for u G (0, 1). Assume that w(u) < for u G (0, 1). 
Then the convexity of R and ^(u,w(u)) > would imply that R(u, 0) > for 
some u G (0, 1). This is impossible by Lemma [231 and we have shown l3~3"b . Finally 
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13.3b follows directly from the definition of w(u) as w(u) — inf {w : (u, w) G L} and 
I as the effective domain of w(u). □ 

We are now ready to show our main result on the long-term properties of the 
log-price process (X t ) t > . 

Theorem 3.4. Let (X t ,Vt)t>o be an affine stochastic volatility model and suppose 
that x(0) < and x(l) < 0. Let w(u) be given by Lemma 3.2 and define 

h(u) = F(u, w(u)), J — {u£ I : F(u, w(u)) < oo} . 

Then [0, 1] C J C I; w(u) and h(u) are cumulant generating functions of infinitely 
divisible random variables and 



lim vj}(t, u, 0) = w(u) for all u G I ; 

t — >oo 

lim — (f)(t, u, 0) = h(u) for all u G J . 

t— >oo t 



(3.5a) 
(3.5b) 

Corollary 3.5. Under the conditions of Theorem \3.4\ the following holds: 
(3.6a) 

U£[0,1 

' 1 



sup \ip(t,u,0)-w(u)\ <Cexp(-X-T) 

!»£[0,1] 



(3.6b) 



sup 

ue[o,i] 



<p(t,u,0) - h(u) 



for some constant C , and with 



< f2Cexp(-£ • T) ; 



d 



X= inf |x(u)| and fl = sup ——F(u,w) 

«e[o,i] ne[0 ,i] ow 



w=0 



Proof. Let u G [0, 1]. By Lenrna|l(ii,Mi(u)) G [0, 1] x (-oo,0]. By Theorem 
F(0,0) = F(1,0) = 0, such that Lemma I2T21 guarantees that [0,1] x (-oo,0] C 
domF. It follows that [0, 1] C J. Define 

z(t, u) = ip(t, u, 0) — w(u) . 

Inserting into the generalized Riccati equation 12. 3bl 

d 

— z(t, u) = R(u, ip(t, u, 0)) = R(u, ij){t, u, 0)) — R(u, w(u)), and z(0, u) = w(u) . 

If if){t, u, 0) < we can bound the right hand side by 

OR 

R(u, ip(t, u, 0)) — R(u, w(u)) < z(t, u)-g—(u, 0) = z(t, u)x(u), 

using convexity of R. By Gronwall's inequality 

z(t,u) < \w(u)\exp(x(u)t) . 
Since x is convex, x(0) < and x(l) < 0, we have shown (13.6a[) . The estimate 



\4>(t,u)-h(u)\ 



(F(u,<f>(s,u)) - F(u,w(u))) ds 



< 



OF 
dw 



(«,0) 



\ip{t,u) - w{u)\ 



yields (|3.6b[) and we have shown Corollary [331 

Let now u G 1° \ [0,1]. Combining Lemma 12.91 and Lemma 13.31 we have that 
R(u,w) > for all w £ [Q,w(u)), and R(u,w(u)) = 0. It follows that the initial 
value ip(Q, u, 0) = is in the basin of attraction of the stable equilibrium point w(u) 
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and thus that ip(t, u, 0) is strictly increasing and converging to w(u). An additional 
argument may be needed at the boundary of /: Let u+ = sup / and assume that 
u + G / (i.e. / is right-closed). Since (u + ,w) G domi? for all w < w(u+), we 
can define §^(u + ,w(u + )) at least as a limit for w | w(u + ). By Lemma \3 . 2 1 either 
w(u + )) < or §p^(u + , w(u+)) — 0. In the first case we can argue as in the 
interior of / that w(u+) is an asymptotically stable equilibrium point. In the second 
case we use once more that by Lemma 12.21 R(u + , w) is, as a function of w, either 
strictly convex or affine. If it is affine, it must be equal to 0, and thus i?(it+, 0) = 0, 
in contradiction to Lemma l2.9l Hence it is strictly convex, and attains its minimum 
at w(u+). This implies that R(u+,w) > for all w G [0, w(u+)) and we conclude 
that ip(t, tt+,0) converges to For u_ = inf J, a symmetrical argument 

applies. 

Assertion (|3.5bj) follows immediately from the representation (12. 4p . and 

1 If 1 

lim —cj)(t,u,0)= lim - / F(u,ip(s,u,0)) ds — F(u,w(u)) 

t^OO t t— »00 t Jq 

for all u G J. 

We have shown that the sequence of infinitely divisible cumulant generating func- 
tions ip(t,u,0) converges on / to a function w(u) that is continuous in a right 
neighborhood of 0. This is sufficient to imply that w(u) is agai n the cumula nt gen- 



erating function of an infinitely divisible ra ndom varia ble (See iFellerl [1971L VIII. 1, 



Example (e)] for the convergence part, and ISato 1999, Lemma 7.8] for the infinite 



divisibility.). The same argument can be applied to (f> and h(u), and we have shown 
Theorem [311 □ 

4. Moment explosions 

In this section we continue to study the time evolution of moments E[5 t u ] = 
E[e x * u ] of the price process in an affine stochastic volatility model. We are in- 
terested in the phenomenon that in a stochastic volatility model, moments of the 
price process can explode (become infinite) in finite time. For stochastic volatility 
models of the CEV-type - a class including the He ston model, but no models with 
jum ps - moment explosions have been studied by lAndersen and PiterbaTel [20071 ] 



an dlLions and Musielal [20071 ]. In the context of option pricing, an interesting result 
of lLed [2004f connects the existence of moments of the stock price process to the 



steepness of the smile for deep in-the-money or out-of-the-money options. Our first 
result shows that in an affine stochastic volatility model a simple explicit expression 
for the time of moment explosion can be given: 

4.1. Moment explosions. By definition, the u-th moment of S t , i.e. E^"] is given 
by Sq exp (</>(t, u, 0) + Vo ?/>(£, u, 0)). We define the time of moment explosion for 
the moment of order u by 

T*(«) = sup{t : E^ 11 ] < oo} . 

It is obvious from the Markov property that ~&[Sf] is finite for all t < T*(u) and 
infinite for all t > T*(u). As in the previous section, the main result follows from a 
qualitative analysis of the generalized Riccati equations (|2.3p . 



^Even though ip(t, u+, 0) converges to note that w(u+) is not a stable equilibrium point in 

the usual sense. This is due to the fact that solutions from a rig/ii-neighborhood N n (io(u+), oo) 
will diverge from w(u+) to +oo. 
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Theorem 4.1. Suppose the conditions of Theorem \3.1\ ho\d. Define J = {u G I : F(u, w(u)) < oo}, 
f+(u) :— sup {w > : F(u, w) < oo} , 
r+(u) := sup {w > : R(u,w) < oo} , 

and suppose that F(u,0) < oo, R(u,0) < oo and x( u ) < °°- 

(a) If u £ J , then 

T*(u) = +oo . 

(b) IfueM\J, then 

rrain(/+(u),r+(u)) ^ 



T.(t») 



i?(w, 77) 



// 0) = 00, i?(it, 0) = 00 or = 00 ^ e?l 

= . 

Proof. Suppose that u G J. Then Theorem 13.41 implies that both ^(t, u, 0) and 
</>(£, u, 0) are finite for all t > 0. This proves (jlj). Let now u G M \ J, 0) < 00, 
i?(w, 0) < 00 and < 00 ■ To prove Jb| we start by analyzing the maximal 

lifetime of solutions to the generalized Riccati equation 

(4.1) J^(t, u, 0) = R(u, i>{t, u, 0)), ^(0, u, 0) = . 

Define M = [0,r+(u)) and note that .) G C(M). Since u £ [0, 1], Lemma l2~9l 
implies that i?(u,0) > 0. It is clear, that at least a local solution tp(t,u,0) to the 
ODE exists, which satisfies < i/j(t, u, 0) < r + (u) and is an increasing fu nction of t 
as lon g as it can be continued. Using a standard extension theorem (e.g. iHartmanl 
[l982l Lem. 1.3.1]) the local solution ip(t, u, 0) has a maximal extension to an interval 
[0,T(u)), such that one of the following holds: 

(i) T(u) =00, or 

(ii) T{u) < 00 and ip(t,u,0) comes arbitrarily close to the boundary of M, i.e. 

lim sup ip(t, u, 0) = r + (u) . 

t->T(u) 

Consider case Since ip is increasing, its limit for t — > 00 exists, but can be 
infinite. Suppose limj^oo ip(t) = a < 00. Then a must be a stationary point, i.e. 
R{u, a) = 0, but this is impossible by Lemma 13.31 The case that a — 00 is only 
possible if r + (u) — 00, such that in this case \im t ^T(u) ipif, u, 0) = r+(u). Consider 
case ([iljl. Since tp is increasing the limes superior can be replaced by an ordinary 
limit and we get limt_>T(«) "4>(t, u, 0) = r + (u) as before. 

Let now T n be a sequence such that T„ "f T(u). By (|2.5p it holds that 
(4-2) / -5T—sds = T n . 

Jo R{u,v) 

R{U,T)) ' 



Letting n — > 00 we obtain that T'(u) = J" r+ ^ ■j^-rds. 
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We can write the time of moment explosion T* (u) as the maximum joint lifetime 
of (j)(t,u,0) and ip(t,u,0), i.e. T*(u) = sup{i > : <fi(t,u,0) < oo A 2p(t,u,0) < oo}. 
By the integral representation (|2.4p it is clear that if f+(u) > r + (u), (f>(t,u,0) 
is finite whenever ip(t,u,0) is finite and T*(it) = T(u). If f+{u) < r+(u) then 

(it), u, 0) = f+(u). Inserting into the representation (14. 2[) yields (|b|. 
For assertion (c), let F(u, 0) = oo, 0) = oo, or x(u) = oo. In the first case, 
<p(t,u,0) does not exist beyond t = 0. In the other cases no local solution to the 
generalized Riccati equation (|4.ip exists, such that ?/>(i,w, 0) explodes immediately. 

□ 



4.2. Moment explosions in the stationary variance regime. In Section 13711 
we have introduced (X t , Vt)t>o as the model in the stationary variance regime. The 
moment explosions of this process can be analyzed in a similar manner as above. 
We define the time of moment explosion in the stationary variance regime by 



Tf (u) := sup |T > : E[S%] < oo j ; 



the superscript 'S' stands for 'stationary'. 

The analogue to Theorem 14. II is the following result: 

Theorem 4.2. Suppose the conditions of Theorem \3.4\ hold. Define f + (u),r + (u) 
as in Theorem\4-.l\ and in addition 



1+ := sup{u> > : l(w) < oo} . 

Suppose that F(u, 0) < oo 7 R(u, 0) < oo and x(0) < oo. 

(a) If u € J and w(u) < Z+, then 

T*(u) = +oo . 

(b) If u £l\ J or w(u) > l +: then 

/■min(/ + (u),r+ («),(+) 

/o ' 

// ^(m, 0) = oo, R(u, 0) = oo or x(0) = oo, then 
(c) 

T?(u) = . 

Corollary 4.3. Under the conditions of Theorem \4-2\ 

T*(u) < T*(u), for allueR 

Proof. By equation (|3.2|) . the moment EfS""] is given by 

E[S?] = exp u, 0) + l(ip(t, u, 0))) . 

This expression is finite, if 4>(t,u,0) and ip(t,u,0) are finite, and if ip(t,u, 0) < Z+. 
It is infinite if 4>{t,u, 0) or ip(t,u, 0) are infinite, or if ip(t,u,0) > Z+. The rest of 
the proof can be carried out as for Theorem 14.11 Note, that now even for u G J, 
the moment can explode, if 1+ is reached by ip(t,u,0) before the stationary point 
w(u). Corollary 14.31 follows easily by comparing the range of integration and the 
conditions for case (a) and (b) between Theorem 14.11 and Theorem 14.21 □ 
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5. Applications 

5.1. Smile behavior at extreme strikes. In the preceding section, we have kept 
u fixed, and looked at the first time T*(u) that the moment E[5"] becomes infinite. 
It will now be more convenient to reverse the roles of T and u, and for a given time 
t to define the upper critical moment by 

u+{t) = sup{u > 1 : E[SJ*] < oo} = sup{u > 1 : T*(u) < t} , 

and the lower critical moment by 

u-(t) = inf {u < : E[Sf] < oo} = inf {u < : T*(u) < t} . 

It is seen that u_ (T) and u+ (T) can be defined as the generalized inverse of T* (u) 
on (—oo,0] and [1, oo) respectively. In addition it is easily derived from Jensen's 
inequality, that 

E[5"] <oo for all u € (u-(t), «+(*)), and 
E[5 t "] =oo for all u € R\ [u_(i), u+(t)] . 
The results of Eel [2004j relate the explosion of moments to the 'wing behavior' 



of the implied volatility smile, i.e. the shape of the smile for strikes that are deep 
in-the-money or out-of-the-money. To give a precise statement, let £ be the log- 
moneyness, which for a European option with time-to-maturity T and strike K is 

given by £ = log(^£^)- 

Proposition 5.1 (Lee's moment formula). LetV(T, £) be the implied Black- Scholes- 
Variance of a European call with time-to-maturity T and log-moneyness £. Then 

V(T,0 _ ?(-u_(T)) 



and 



limsup 

y(T,Q ?MT) - 1) 
limsup — — = 



where c(x) = 2 — 4 (v 'x 2 + x — x) and u±(T) are the critical moment functions. 

The function q is strictly decreasing on R^oi mapping to 2, and oo to 0. 
Thus for fixed time-to-maturity T, the steepness of the smile is decreasing with 
|w±(T)|. A finite critical moment u±(T) implies asymptotically linear behavior of 
V(T, £) in ^, and an infinite critical moment implies sublinear behavior of V(T, £). 
It is also evident that U-(T) determines the 'left' side of the volatility smile, also 
known as small-strike, in-the-money-call or out-of-the- money-put side; u+(T) de- 
termines the 'right' side, or large-strike, out-of-the-money-call, in-the-money-put 
side. Finally we mentio n that Lee's result has been extended and strengthened by 



Benaim and Fria [2006] from a 'limsup' to a genuine limit under conditions related 



to regular variation of the underlying distribution function. 

5.2. Forward-smile behavior. The forward smile is derived from the prices of 
forward-start options. For a forward-start call option - all options we consider 
are European - a start date r, a strike date T + r and a moneyness ratio M are 
agreed upon today (at time t = 0). The option then yields at time T + r a payoff 

of — A'l^j , i.e. the relative return over the time period from t to r + T, 
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reduced by M and floored at 0. Under the pricing measure the value of such an 
option at t = is given by 



(5.1) 



-r(T+r) E 





= e- Tr E 







where we define the log-moneyness £ of a forward-start option as £ = log M + rT. 
Forward-start options are not just interesting in their own right, but ar e used as 
buildi ng blocks of more complex derivatives, such as Cliquet options (see iGatherall 
[20061 Chapter 10]). 

Analogously to plain vanilla options, we can define the implied forward volatility 

cr(r, T, £), by comparing the forward option price to the price of an option with 
identical payoff in the Black-Scholes model. Note that the implied forward volatility 
depends also on r, the starting time of the contract. For r = 0, the implied volatility 
of a plain vanilla option is retrieved. More interesting is the behavior for r > 0. 
Intuitively, we expect the implied volatility (and the option price) to increase with 
t in a stochastic volatility model, since the uncertainty of the variance V T at the 
starting date of the option has to be priced in. In an affine stochastic volatility 
model, it will be seen that under mild conditions, the implied forward volatilities 
<t(t, T, £) actually converge to a limit as r — > oo. Not surprisingly, this behavior 
is related to the convergence of (Vt)t>o to its invariant distribution. In the limit 
t — ► oo , the pricing of a forward-start option is equivalent to the pricing of a plain 
vanilla option in the stationary variance regime (cf. Section 13. lj) . 



Proposition 5.2. Let (Xt, Vt)t>o be an affine stochastic volatility model, satisfying 
the conditions of Proposition HOI Let <t(t, T, £) be the implied forward volatility in 
this model. Then 

km a(r,T,0 = a(T,0 , 



= A T 



where a(T, £) is the implied volatility of a European call with payoff ^< 
and Xt is the log-price process of the model in the stationary variance regime. 
Proof. We can write the price of a forward-start call as 



C(r,T,0 = e- Tr E 



Xt+t — X t 



E 



Denote by C (T, £, er) the (plain vanilla) call price in a Black-Scholes model with 
volatility a and the normalization So = 1. It is easy to see that the price of a 
forward-start option in the Black-Scholes model is just the discounted plain vanilla 
price, i.e. C bs (t,T : £,ct) = e~ rT C BS (T,$,,cr). By definition, the implied forward 
volatility of the call C(r, T, £) satisfies 



C BS (T, & <t(t, T, 0) = e rT C(r, T, £) = E 
Taking the limit r — * oo on both sides we obtain 



E (0,K) 



C BS (T,£, lim a(r,T,£))=E 



lim E (0 ' K 



(e XT 



= E 



e x - - e< 



using dominated convergence. It is well known that the above equation allows a 
unique solution in terms of the Black-Scholes implied volatility, and we get <j(T, £) = 
lim T _ KX) (t(t,T,Z). □ 
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Combining Lee's moment formula with our results on moment explosions under 
the stationary variance regime (Theorem 14. 2j) . asymptotics of a(T,£) for £ — > ±00 
can be derived. 

6. Examples 

6.1. T he Heston model with and without jumps. In the model of iHestonl 



1993], the log-price (X t )t>o and the corresponding variance process (V t )t>o are 



given under the risk-neutral measure by the SDE 



dX t = ~dt + y/VidW? 
dV t = —X(Vt -6)dt + C^/v t dW t 



where , are Brownian motions with correlation parameter p, and (,,\,6 > 0. 
In affine form, the model is written as 

(6.1a) F(u,w)=\6w 

1 c 2 

(6.1b) R(u, w) = — (u 2 — u) + ~2 w2 — ^ w + UW P( ■ 

It is easily calculated that \ is given by xi u ) = PC U ~ A. We will first analyze 
the long term behavior of (X t ) t > , with the help of Theorem 13.41 To satisfy the 
condition x(l) < we need A > (p. Note that this condition is always satisfied if 
p < 0, the case that is typical for applications. Solving a quadratic equation we 
find that 



w(uj — 1 — , and li(u) = Xuw{u) 



(A - upC) - y/(\-up( f-e(u2-u) 

c 

Denoting the term under the square root by A(u), we see that w(u) and h(u) are 
both defined on J = I = {u : A(it) > 0}. Since R is a second order polynomial in 
the Heston model, the equilibrium points of the generalized Riccati equation for ip 
form an ellipse in the (u, w)-plane, and w(u) is given by its lower part - see Fig- 
urc[T]for an illustration. Interestingly, w(u), and also h{u), are cumulant generating 
functions of a Normal Inverse Gaussian distribution (cf. iBarndorff-Nielseii 1997 



Eq. (2.4)]). Thus, for large t, the price process of the Heston model is, in terms 
of its marginal distributions, close to a Normal-Inverse-Gaussian exponential-Levy 
model. 

Next we consider moment explosions in the Heston model. As mentioned above, 
moment ex plosions in the Heston model (and other models) have already been 
studied by lAndersen and Piterbargl [20071 ] . Nevertheless this will provide a first 
test of Theorem 14.11 In the case of the Heston model it is easily determined from 
(16. ip that f+{u) = r + (u) = 00. Calculating the integral in case © of Theorem HU 
we obtain 

( +00 A(u) > 

(6.2) T,(t*) = <^ 2 / , , \ . , , . 

1 -^rm ( v arctan + 7ri ix(»)<o} j a (u) < . 

In Figure [Ha plot of this funct ion for typical parameter values is shown. Note that 
lAndersen and Piterbargl 2007 distinguish an additional case where x( u ) < 0> Du t 



A(u) > 0. A little calculation shows that this can only happen if x(l) > 0, a case 
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u 

FIGURE 1. This plot shows the stable and unstable equilibria of the gen- 
eralized Riccati equation of a Heston model with parameters p = —0.7165, 
C = 0.3877, A = 1.3253 and 6 = 0.0354 (taken from iGatherall [20061 Ta- 
ble 3.2]). It can be seen how the solutions ip(t,u) converge to the stable 
equilibrium points, which form the lower boundary of an ellipse in the (u, w)- 
plane. 



that is precluded by our assumptions in Theorem l3.4i and never occurs when p < 0. 

We will now study the effect of adding jumps to the Heston model. The simplest 
case is the addition of an independent jump component with constant activity: 
Let (Jt)t>o be a pure-jump Levy process, independent of (M^ t 1,2 )t>o and define the 
Heston-with-jumps model by 

dX t = fs - dt + y/VidW? + dJ t 

dV t = -X(V t -9)dt + (y/VtdW? . 

The drift S is determined by the martingale condition for (S t )t>o- To make the 
example simpler, we assume that ( Jt)t>o jumps only downwards. This is equivalent 
to saying that the Levy measure m(dx) of (Jt)t>a is supported on (— oo,0). The 
affine form of the model is 

(6.3) F(u,w) = \6w + k{u) 

1 c 2 

(6.4) R(u, w) = — (u 2 — u) + —w 2 — Am + uwpQ , 
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where k(u) is the compensated cumulant generating function of the jump part, i.e. 



k(u) = / (e xu - 1) m(dx) - u I {e x - 1) m{dx) . 

Let K- < be the number such that k(u) is finite on (k_, oo) and infinite outside. 
For example, if the absolute jump heights are exponentially distributed with an 
expected jump size of 1/a, then k_ = —a. 

To analyze the explosion times of this model, note that R, and thus x{ u )i w i u )i I 
and r + (u) have not changed compared to the Heston model. As long as u > 
the explosion time T*(u) is the same as in the Heston model. However, if u < k_, 
F(u,0) = oo and by Theorem 14. li T*(it) = 0. Thus, the addition of jumps to 
the Heston model has the effect of truncating the explosion time to zero, whenever 
u < K-. 

From the viewpoint of the critical moment functions, does not change com- 

pared to the Heston model, but u—(t) does; in the model with jumps it is given 

by 

U i ump (t) = u Hcston (i) Vk_ . 
Since u_ is increasing with f, it makes sense to define a cutoff time Xjj by 
T t = sup {t > : u Hoston (i) =K-} = T,(k_) , 

such that 

w Hcston (t) < ui ump (t), if t < Tj 
u Hoston (<) = ui ump (t), if t > T t . 

In Figure [2] a comparison of the critical moment functions in the Heston model with 
and without jumps can be seen. By Lee's moment formula, the critical moment 
u_(t) moving closer to will cause the left side of the implied volatility smile to 
become steeper. Thus the net effect of adding the jump component (Jt)t>a to the 
Heston model, is a steepening of the left side of the smile for maturities smaller 
than T(j. For times larger than Xjj, the asymptotic behavior of the smile (in the 
sense of Lee's formula) is exactly the same as in the Heston model without jumps . 



sense oi .lcc s iormuiaj is exactly tnc same as m tne rieston model witnout jumps. 
This corresponds well to the frequently made observation (see e.g. iGatheral 20061 



Chapter 5]) that a Heston model with jumps can be fitted well by first fitting a 
(jump-free) Heston model to long maturities, and then calibrating only the ad- 



ditional parameters to the full smile. In fact iGatherall proposes (on heuristical 
grounds) the concept of a 'critical time' T, after which the influence of an inde- 
pendent jump component on the implied volatility smile can be neglected. The 
analysis of the Heston model with jumps is of course easily extended to the case 
that (Jt)t>o is not one-sided. In that case the effects discussed above will be seen 
to affect also the right side of the implied volatility smile. 

6.2. A model of Bates. We consider now the model given by 

dX t = (s - dt + \/VtdWt + J x N(V t ,dt, dx) 



dV t = -X(V t -0)dt + (VVidWt 



where as before A, 9, £ > and the Brownian motions are correlated with correlation 
p. The jump component is given by N(Vt,dt,dx) = N(Vt,dt,dx) — n(Vt,dt,dx), 
where N(Vt, dt, dx) is a Poisson random measure, and its compensator n(Vt, dt, dx) 
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is of the state- dependent form Vtp(dx)dt, with u(dx) the Le vy measure given in 
(|2.6|) . A model of this kind has been proposed by iBatesI 2000| to e xplain the time- 



variation of jump-risk implicit in observed option prices. Bates! also proposes a 
second variance factor, which we omit in this example, in order to remain in the 
scope of Definition ^. 81 It would however not be difficult to extend our approach to 
the two-factor Bates model, since the two proposed variance-factors are mutually 
independent, causing the corresponding generalized Riccati equations to decouple. 
Since it is affine, the above model can be characterized in terms of the functions F 
and R: 

(6.5) F(u,w) = X9w 

1 c 2 

(6.6) R(u, w) = — (u 2 — u) + ~2 w2 ~ ^ w + uwp( + n(u) . 

where k(u) is the compensated cumulant generating function of the Levy measure 
p. As in the Heston model we can obtain w(u) and h(u) explicitly, and get 



h(u) — — ^ - ^ - - , and h(u) = X9w(u) , 
C 

where x( u ) — p( u ~ A and A(u) = x( u ) 2 ~ C 2 ( m2 — u + 2/c(u)). Both w(u) and h(u) 
are defined on / = J = {u : A(u) > 0}. The time of moment explosion can again 
be calculated explicitly, and is given by 

!+oo A(it) > 

__J_ ^ arctan + 7rl {x{u)<0} ^ -oo < A(u) < 

A(u) = -oo . 

6.3. The Barndorff-Nielsen-S hephard model. The Barndorf f -Niels en-Shephard 
(BNS) model was introduced bv lBarndorff-Nielsen and Shephardl 2001 1 as a model 
for asset pricing. In SDE form it is given in the risk-neutral case by 

dX t = (S- ^V t )dt + ^VtdWt +pdJ xt 
dV t = -XV t dt + dJxt 

where A > 0, p < and (Jt)t>o is a Levy subordinator, i.e. a pure jump Levy 
process that increases a.s. The drift 5 is determi ned by the martingale condition 
for (S t )t>o- The time-scaling J\ t is introduced bv lBarndorff-Nielsen and Shephardl 
to make the invariant distribution of the variance process independent of A. The 
distinctive features of the BNS model are that the variance process has no diffusion 
component, i.e. moves purely by jumps and that the negative correlation between 
variance and price movements is achieved by simultaneous jumps in (V t ) t >o and 
(X t )t>o- The BNS model is an affine stochastic volatility model, and F and R are 
given by 

(6.8) F(u, w) — Xk(w + pu) — uXk(p) 

(6.9) R(u, w) = i(u 2 - u) - Xw 
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where k(u) is the cumulant generating function of (Jt)t>o- 

We simply have x( u ) = — A and w(u) from Lemma 13.21 is given by 

w(u) = ^-(u 2 -u) . 



It follows that 



//!//) : : A;,- ( ^ +u I p- — ) j //A/vl/.i 



This expression can be interpreted as cumulant generating function of a Brownian 
motion with variance j- and drift p — -j , subordinated by the Levy process J\t and 
then mean-corrected to satisfy the martingale condition. 

To analyze moment explosions in the BNS model, let n+ := sup {u > : k(u) < oo}. 
It is easy to see that / + is given by / + = max(« + — pu, 0). Since r + = oo, we have 
that the explosion time for the moment of order u is given by 

f U Jn = _l log ( 1 2A(max(. + -p U ,0)) 



o R(u, j]) A \ u(u — 1) 

The critical moment functions u±(T) can be obtained explicitly by solving a qua- 
dratic equation, and are given by 




The large-strike asymptotics for the implied volatility smile in the sense of Lee can 
be explicitly calculated by inserting u± into Proposition ^. II 



6.4. The Heston model in the stationary variance regime. In the Heston 
model the limit distribution of the variance process (Vt)t>o is a Gamma distribu- 
tion with parameters (— ^r, -p-). This is well-known, but can also be obtained by 
applying Proposition [3Tj The cumulant generating function l(w) is thus given by 



l(w) 



2X9 
"7 



, log ( 1 



c 2 

2T 



defined on (— oo, #), such that I 



,,^u i + — || . As before we have that x( u ) — pC u ~ A, and 
we assume that x(l) < 0. In addition we define x + ( u ) = p(u + \. By Theorem 14. 21 
the explosion time in the stationary regime is given by 



(6.10) T?(u) 



2A/C 2 



-A 



< + X +2Ay / A-A 



C+X-2AVA-A 



VAw)>-x + W, 



0< ^a)<- X + (u), 



arctan + 7rl {x+x<A} ) A(«) < . 



In Figure [5] Tf(u) is plotted together with T*(u) for the Heston model. 
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FIGURE 2. This plot shows the critical moment functions u±(t) for a Heston 
model with the same parameters as in Figure [l] Also shown are u? (t) for 
the model in the stationary variance regime, and u^ ump (i) for the Heston 
model with an independent jump component, whose negative jump heights are 
exponentially distributed with mean a = —0.1. Note that u^_ ump (t) coincides 
with u±(t) everywhere except in the lower left corner of the plot. 



6.5. The BNS model in the stationary variance regime. In the BNS model, 
the cumulant generating function of the limit distribution L of the variance process 
is given by Proposition 13.11 by 

l(w) = / d-q , 

Jo V 

provided the log-moment condition J y>1 (logy) fi(dy) < oo holds for the Levy mea- 
sure of (Jt)t>o- The above integral is finite as long as w G (— oo, and infinite 
outside. Thus /+ = «+. In Section 16.31 we obtained that f+(u) = k+ — pu, such 
that the time of moment explosion under stationary variance is given by 

^ , f< 2Xk{u) 



R(u, rj) A \ u(u — 1) 

where k(u) = k+ for u > 1 and k(u) = max(«; + — pu, 0) for u < 0. Again, this 
expression can be inverted to give the critical moment functions in the stationary 
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variance case. By definition p < 0, such that we obtain 



u b .{T) 



4(T) = 



1 - e 



-AT 



1 (2k+ - p)X 
4 1 



p 2 X 2 



-AT 



(1 



,-AT 



2k , A 



2 ' V 4 1 - e" AT ' 
Appendix A. Additional proofs 
Proof of Theorem \2.1\ Let t < r. By the flow equation we can write 
</>(t, u, 77) = <£(i, it, 77) + <Xt - i, u, ^(i, u, 77)) 
ip(r, u, 77) = ip(r - t, u, ip(t, u, 77)) . 

Since the left sides are finite by assumption, it follows that also <p(t, u, 77) and 
i/j(t,u,r}) are. Vt is non-negative, such that 

|E [cxp (uX t + wV t )}\ < |E [exp (uX t + r,V t )]\ , 

whenever Rew < Re 77. Thus (j>(t,u,w) and %l)(t,u,w) exist for all w € C with 
ReTjj < Re 77. As a particular case we can conclude that (f>(t,u,w) and ip(t,u,w) 
exist for all (u,w) in U :— {(u, w) £ C 2 : Re 77 = 0, Rew < 0}. 

We also define U° := { (u, w) e C 2 : Re u — 0, Re w < 0} , and show next that <f)(t, 77, ttj) 
and ip(t,u,w) are (right-)difier entiable at rj = for all (u w) 6 The key idea 
of our proof is originally due to [Montgomery and Zippinl fl955t. and ha s also been 
presented in iFilipo vie and Teichmannl 120031 and Daws on and Lil First note 

that the identity 

d d \ 

— <f)(t, u, w) + Vo—ipit, u, w) exp (<j)(t, u, w) + V ^(t, u, w) + 



E \wV,e uXt+wVt ] = 



dw 



dw 



shows that -^<p(t, u, w) and ^t/j(t;, u, w) exist, and are continuous for all t < t and 
(u, w) 6W°. By Taylor expansion it holds that 



ip(r, 7i, 7/1(7;, u, w)) dr - 



1 d 

tfj(r,u,w) dr — I —— ip(r,u,w) dr (ij}(t,u,w) — w) 
ow 



(A.l) + o(\ip(t,u,w) - w\) 

On the other side, using the flow property, we calculate 



tp(r, u, 7/1(7;, u, w)) dr - 



ijj(r, 77, ttj) dr 



%p(r + t, 77, ttj) dr - 



ip(r, 77, ttj) dr 



(A.2) 



s+t 



ip{r, 77, ttj) dr - 



vp(r,u,w)dr= / ip(r + s, u, w) dr ■ 



ip(r, u, w) dr . 



Denoting the last expression by I(s,t), and putting (|A.f [) and (|A.2|) together, we 
obtain 



lim ■ 



d 



*->0 \lp(t, 77, 7«) — 77j| 

Thus, writing M s — ^ J Q S -J=^ip(t, u, w) dr, we have 



—ip(t,u,w)dr 



lim - \tp(t, 77, ttj) — w\ 



I(s,t) 



lim 

t^o st 



■ \M S 



ip{s, u, w) — w 



\M S 
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But M s is a continuous function of s, and lim s _^o Ms = ^^(O; M i w ) = 1; such that 
for s small enough M s ^ 0. We conclude that the left hand side is finite, and using 
(jA.ip we obtain that 



lim 

t->o 



tp(t,u,w) — W ( ip(s,u,w) — w 



t 



d_ 

dw 



ip(r, u, w) dr 



The finiteness of the right hand side implies the existence of the limit on the left. 
In addition the right hand side is continuous for (u, w) £ U°, showing that also the 
left hand side is. A similar calculation for tfi(t, u, w) shows that 



lim 

t->o 



(p(t,u,w) 4>{s,u,w) 



( ip(t, u,w) — w 

lim 

t->o \ t 



d_ 

dw 



(r, u, w 



dr 



allowing the same conclusions for cf>(t,u,w). We have thus shown that the time- 
derivatives of <j)(t, u, w) and ib (t, u, w) at t = exis t, and are continuo us in U° . 
Combining Duffie et al. 20031 Proposition 7.2] and Du ffie et al. 2003 . Proposi- 
tion 6.4] the differentiability can be extended from U° to U, and we have shown 
th at (Xt, Vt)t> o is a regular affine process. The rest of Theorem 12.11 follows now as 
in lDuffie et all [20031 Theorem 2.7] □ 



Proof of Lemma \2.SX We prove the assertions of Lemma 12.21 for F; they follow 
analogously for R. By the Levy-Khintchine representation (|2.6|) . F(u, w) + c is the 
cumulant generating functions of some infinitely divisible random variables, say X. 
Writing z — (u,w) € K 2 , and using Holder's inequality it holds for any A e [0,1] 
that 



(A.3) F(\ Zl + (1 - X)z 2 ) = logE 



e \(z lt X) e (l-\)(z 2 ,X) 



C < 



< A log! 



(1 - A)E 



c = XF( Zl ) + (1 - X)F{z 2 ) , 



showing convexity of F. In addition equality in (|A.3|) holds if and only if ke^ Zl ' X ^ — 
e {z2,x) a s f or some k > 0. This in turn is equivalent to (zi — z%, X) being constant 
a.s. Choosing now z\ and z 2 ^ z\ from some one-dimensional affine subspace 
U = {p+ (q,x) :iel] of R 2 , we see that either (q,X) is constant a.s. in which 
case F\u is affine, or it is not constant, in which case strict inequality holds in (|A.3|) 
for all zi,Z2 € U, showing (c). 

Let L a — {z : F{z) < a} be a level set of F, and z n € L a a sequence converging to 
z. Then by Fatou's Lemma 



logE[e <z ' X> ] - c < liminf logE[e 



■ A') 



— c < a 



showing that z £ L a and thus that F is a closed convex function. Finally F is 
proper, because F(0,0) = c > — oo, showing (a). 

Next we show analyticity: Consider the random variables X n :— Xlnx\< n y- Since 
they are bounded, their Laplace transforms, and hence also their cumulant gener- 
ating functions are entire functions on C 2 , and thus analytic on R 2 . As a uniform 
limit of analytic functions F(u, w) is analytic in the interior of domF, showing (b). 
Assertion (d) follows directly from Theorem 12. 11 □ 
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